Seeing the orbital ordering in Iron-based superconductors with magnetic anisotropy 
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The orbital fluctuation of the conduction electrons in the Iron-based superconductors is found to 
contribute significantly to the magnetic response of the system. With the use of a realistic five-band 
model and group theoretical analysis, we have determined the orbital magnetic susceptibility in 
such a multi-orbital system. At n = 6.1, the in-plane orbital magnetic susceptibility is predicted 
to be about Wfj,%/eV, which is more than 2/3 of the observed total susceptibility around 200 
K in 122 systems(of about 14/ig/eV or 4.5 x 10 -4 erg/G 2 moiA3— ). We find the in-plane orbital 
magnetic response is sensitive to the breaking of the tetragonal symmetry in the orbital space. In 
particular, when the observed band splitting(between the 3d xz and the 3d HZ -dominated band) is 
used to estimate the strength of the symmetry breaking perturbation^, a 4.5% modulation in the 
in-plane orbital magnetic susceptibility can be produced, making the latter a useful probe of the 
orbital ordering in such a multi-orbital system. As a by product, the theory also explains the large 
anisotropy between the in-plane and the out-of-plane magnetic response observed universally in 
susceptibility and NMR measurements. 

PACS numbers: 



An unresolved issue in the study of the Iron-based su- 
perconductors is the role of their multi-orbital nature. In 
most other superconductors, the orbital degree of free- 
dom is quenched at low energy in the crystal field envi- 
ronment. However, both LDA calculation and ARPES 
measurements'— indicate that in the Iron-based super- 
conductors all the five Fe 3d orbital play essential role 
in forming the low energy degree of freedom around the 
Fermi surface. Many novel properties of the Iron-based 
superconductors, especially those in the name of elec- 
tronic nematicity^ , have been argued to be related to 
the orbital ordering in these systems*^—. Most recently, 
a two-fold modulation of the magnetic susceptibility in 
the Fe-Fe plane is found to develop around a tempera- 
ture that is significantly higher than the structural phase 
transition point 5 . However, it is still a mystery how the 
observed electronic nematicity is related to the orbital 
ordering of the system. 

Another puzzle about the Iron-based superconduc- 
tors is the strong anisotropy in their magnetic response 
observed universally in susceptibility and Knight shift 
measurements*^—. The susceptibility in the Fe-Fe plane 
is found to be significantly larger than that perpendicular 
to it. This is very unusual, since the magnetic response 
of a transition metal is usually attributed to the spin of 
its valence electron and is essentially isotropic. The or- 
bital magnetic response, on the other hand, is usually 
quenched as a result of the crystal field effect. However, 
since the crystal field splitting in the Iron-based super- 
conductors is very small and all the five 3d orbital are 
involved in the low energy physics^—, the orbital angu- 
lar momentum of the conduction electron can contribute 
to the magnetic response of these systems. Such a con- 
tribution is intrinsically anisotropic and depends on the 
electronic structure of the system, especially on the sym- 
metry breaking in the orbital space. 

The purpose of this paper is to evaluate orbital mag- 



netic response of the Iron-based superconductors from 
a realistic model and to explore the relation between 
orbital ordering and the electronic nematicity observed 
in recent torque magnetometry measurement^. We find 
the orbital magnetic susceptibility in these multi-orbital 
systems is comparable in magnitude with the measured 
total magnetic susceptibility. More specifically, the in- 
plane orbital magnetic susceptibility is predicted to be 
about 10/ig/eV, which accounts for more than 2/3 of the 
observed susceptibility at 200 K in 122 systems* ^ 16 ' 22 . 
Furthermore, the in-plane orbital magnetic response is 
found to be sensitive to the breaking of the tetragonal 
symmetry in the orbital space, making it a useful probe 
of orbital ordering in these multi-orbital systems. As a by 
product, the observed strong anisotropy between the in- 
plane and out-of-plane magnetic susceptibility also find 
a natural explanation from our calculation. 

The Iron-based superconductors have a very compli- 
cated band structure. In this study, we adopt the five- 
band tight-binding model derived from fitting the LDA 
band structure^ of the LaFeAsO system. Following the 
notations of Ref 20, the band model reads, 
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where v, v' = 1 , . . , 5 is the index for the five maximally 
localized Wannier functions (MLWFs) on the Fe site, 
namely, |1) = \3ekz*-K>), |2) = \3d xz ), |3) = \3d YZ ), 
1 4) = \3d x z-Y 2 ) and |5) = \3dxy)- t^'J denotes the 
hopping integral between the v-th and v'-th orbital at 
site i and site j. Here an unfolded scheme is adopted 
as m Refill The X and Y-axis for the Wannier func- 
tions, which are in the Fe-As bond direction, are rotated 
by 45 degree from the x and y-axis of the Fe-Fe square 
lattice(see Fig [J). is the on-site energy of the z/-th 
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orbital. The hopping integral is truncated at the fifth 
neighbor and the values of the model parameters can be 
found in Refill 




FIG. 1: The square lattice of the Fe ions (shown as gray 
dots) and the local coordinate system for the atomic orbital. 
The red and blue dots denote the As ions above and below 
the Fe-Fe plane. 4> ls the angle between the X-axis and the 
direction in which magnetic susceptibility is measured. In 
the tetragonal phase, the point group symmetry around the 
Fe ion is D2d, which is broken down to D2 in the orthogonal 
phase. 



The interaction of electron has the following general 
form 



Hint = U n i:Vja ni t u,a + (V - J) »>.,:«» 



-U' 2J n i^cr n i 

-J V c f 



(2) 



Here we have included the intra- and inter-orbital 
Coulomb repulsion, the Hund's rule coupling and the 
pair hopping term and have assumed that U 1 = U — 2 J. 



is the number density 



operator of the electron and a = — a. 

Since the five Fe 3d orbital \v) are all real functions, 
they can not carry current and thus their orbital angular 
momentum are quenched in the static limit. However, 
since the orbital content varies on the Fermi surface, 
fluctuation in the orbital character and orbital angular 
momentum survives in the low energy limit and can con- 
tribute to the magnetic response of the system. In the 
following we will calculate such a magnetic response in 
the RPA scheme. 

The orbital magnetic susceptibility is defined through 
the correlation function of the orbital magnetic moment 



in the following way 

X £(q,T)=-(T T L Q (q,T)L"(-q,0)), 



(3) 



in which L a (q, r) denotes the Fourier component of the 
orbital magnetic moment density in the a direction and 
a = X, y, Z . Here we use /x^ as the unit of susceptibility. 
The operator for the orbital magnetic moment on a given 



site is defined as L a = J2uy,a c l,J», 



V,<T V,V' ^ v if' 



where if*. 



is the matrix element of the orbital magnetic moment in 
the basis spanned by the five MLWFs. 

The matrix element i"„, can be determined in prin- 
ciple from a first principle calculation. Here we will be 
satisfied with the result of a semi-quantitative analysis, 
for which much simplification can be achieved when sym- 
metry arguments are adopted. In the following, we will 
illustrate the steps for l z v i- First, since L z is time re- 
versal odd and the five 3d orbital are all real, l z v , must 

be purely imaginary. Second, since L z is odd under the 
action of the three generators of the £>2<i point group 
around each Fe ion, namely R x (ir), ax and ay—, while 
the five 3c? orbital transform as 
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the only none-zero matrix elements are if 3 = — l z 2 and 
l z 5 = —if 4. Thus L z can be generally written as 

L Z = j^h lC 2,(T C 3,«T + I2cl^c 5 . a ] + h.c, 
a 

in which 71 and 72 are two real numbers. Following the 
same line of reasoning one find that 

L X =i X^ 73C l,f C3 ' CT + T 4C 2, CT C 5,ct + 754,<t C 4,<t] + h.c. 
a 

L Y = -i X^ 73C -L C2 ' CT +^ C l.a C 5,a ~ 754,<t C 4,<t] + h.C, 
a 

with the three real coefficients 73,4.5 left undetermined. 
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To have an estimate of the values of the five coefficients 
7i,..,5) we approximate the five MLWFs v = 1, ..,5, 
with the five Fe 3d orbital in the atomic limit. These 
atomic orbital are related to the spherical harmonics of 
I = 2 in the following ways (apart from the radial part of 
the wave function which is not used in determining the 
matrix element of L a ) 



|1) 
|2) 

|3> 

|4> 

|5) 



|2,0) 

7=2 (l 
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V2 
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V2 



(|2,-1) + |2,1)) 
(|2,-2} + |2,2» 



_(|2,-2}-|2,2)), 



where |2,m) oc Y™ are the spherical harmonics oil = 2. 
Since (2, m'\L z \2, m) = mh m ^ and (2, m\L + \2, m') = 
\J& — m'(m' + l)5 mim / + i(here L + — L x +iL Y ), we have 

71 = - 74 = 75 = -1 

72 = -2 

73 = \/3. 

We will use these values in the following calculation. 

The bare orbital magnetic susceptibility is readily ob- 
tained as 



X °£ a (T) = lim I Y 



/(£k+q,m') - / (£k,m| \ ja |2 
| k.m,m' | 



k, m,m' 



Ck+q,-. 
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Here £k,m = £k,m — /U is the band energy of the ru- 
th. band(m = f,...,5) and fi is the chemical potential. 
L Z,rn, m ' = Y, v y l uyK,u,m u ^y,m' and u k ^ m is the wi- 
th eigenvector of the band Hamiltonian at momentum k. 
As a comparison, the Pauli spin susceptibility is given by 
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(T) = lim-V 

q^O N ^ 
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Unlike the orbital magnetic susceptibility, the Pauli spin 
susceptibility has contribution only from intra-band pro- 
cess. Thus at low temperature the spin susceptibility 
is solely determined by the electronic state around the 
Fermi surface, while the orbital magnetic susceptibility 
depends on electronic states both on and far away from 
the Fermi energy. As a result, both the temperature and 
the doping dependence of the orbital magnetic response 
should be much weaker than that of the spin magnetic 
response. 

Now we consider the RPA correction of the orbital 
magnetic susceptibility. The orbital magnetic excita- 
tion of the system has the general form of O vv = 
*£o-(4 a c v' > _ c l>,a c v,<?)- Without losing generality, we 



assume v' > v. There are in total 10 such excitations 
and all of them are time reversal odd and spin singlet. 
The correlation function between these excitations can 
be defined in the following way 



Xo 



(d,r)=-(T T O vv (q,r)0™ (-q,0)), 



and the corresponding bare susceptibility in the static 

limit Xo' 1 ' VV CO 1S gi ven by an expression similar to 
Eq.Q, except that the matrix element |L£ m m ,\ 2 should 
be replaced by 

^ (^k^m^k^'jm' ^k.u'.m^k^m' ) ■ 



The RPA correction of Xo ^ s contributed by the 
inter-orbital Coulomb repulsion, the Hund's rule coupling 
and the pair hopping term. The RPA kernel is extremely 
simple and is given by V vv i }VV / = ^ U 7 J ^ 8 vv > iVV i (see Sup- 
plementary material A). The RPA corrected susceptibil- 
ity can be written formally as 



XO 



Xo 



in which xo, Xo an d V are au to be understood as 10 x 10 
matrix(we note while V is a diagonal matrix in the space 
of O v ' u , Xo 1S n °t)- The orbital magnetic susceptibil- 
ity can be obtained from the combinations of the matrix 
element of xo- For example, 



23,23 

Xo 



+ 4(x 4 o M5 



i 23,45\ 

+ X Q ' )• 



The orbital magnetic susceptibility in other direction can 
be obtained in a similar way. 

The observation of the two-fold modulation in the in- 
plane magnetic susceptibility indicates that the tetrago- 
nal symmetry of the system is broken down to orthog- 
onal. This can happen either through orbital ordering, 
or through nematicity in spin correlation 18,19 . Here we 
assume it happens through orbital ordering, since the or- 
bital magnetic response is much more sensitive to it than 
to spin nematicity. The form of the symmetry break- 
ing perturbation in the orthogonal phase can be largely 
determined by group theoretical arguments. Among the 
five 3d orbital, the 3d 3 z 2 -R 2 , 3dxy and 3d X 2 -y 2 or ~ 
bital each form a one dimensional representation of the 
D2d point group. The 3dxz and 3dyz orbital form a 
two-dimensional representation which becomes reducible 
when the symmetry is lowered to orthogonal. We thus 
focus on symmetry breaking terms in the space spanned 
by the 3dxz and 3dyz orbital. A group theoretical anal- 
ysis then shows that up to nearest neighboring hopping 
terms, the only allowable symmetry breaking perturba- 
tion in the orthogonal phase takes the form (see Supple- 
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mentary information B) 

AH = m^2(cl2,*^,3,<? + C U,a C i,2,a) 
i,cr 

+ m ^2 d s(cl,2^ c t+S,2,a + ct 3)CT C i+ 5 !3!(7 ) 

2,(5,(7 

+ ^3 X^( C l2,<r C ''+<5,3, CT + q,3,<rCi+«,2,<r), (5) 

2,5,(7 

in which 8 = ±x, ±y is the vector between nearest neigh- 
boring Fe sites, ds is the d-wave form factor and d± x = 1, 
d± y = — 1. Here, 771 is the strength of the on-site symme- 
try breaking perturbation. 772 and 773 are the strengths 
of the d-wave intra-orbital and s-wave inter-orbital hop- 
ping terms between nearest neighboring Fe sites. From 
ARPES measurement^, it is found that the splitting be- 
tween the ?>d xz and the 3d yz -dominated band is zero at 
the r point and maximizes at the X and Y point. Among 
the three perturbations in Eq. (|S|) , only the d-wave intra- 
orbital hopping term is consistent with such a momen- 
tum dependence. For example, both the 771 or 773-type 
perturbation would result in an nonzero band splitting 
at the r point, which is not observed. Furthermore, the 
^3-type perturbation has no effect at the X and Y point, 
where the observed band splitting reaches its maximum. 
We thus set 771 =7/3 = 0. This leaves us 772 as the only 
undetermined parameter. 

We are now at the position to present the numerical 
results. Our calculation is done at a fixed band filling 
of n = 6.1. The chemical potential is determined by 
solving the mean field particle number equation at each 
temperature. We have set U = 1.2eV, J = 0.15eV, as is 
chosen in Ref|M To estimate the value of 772 from the 
observed band splitting, we note that the band width 
of the Iron-based superconductors is significantly smaller 
than the prediction of band structure calculation. We 
thus fit the relative rather than the absolute magnitude of 
the band splitting. According to ARPES measurement, 
the maximal band splitting between the id yz and Zd xz - 
dominated band is about one half of the dispersion of the 
Sciy^-dominated band between the T and X poinlji. To fit 
such a splitting, we set 772 = 30 meV. The calculated band 
dispersion along the V — X and T — Y direction is shown 
in Fig|2j which looks very similar to the experimental 
result-. The temperature dependence of 772 is modeled 
by the mean field form of 772 (T) = 772(0)^1 - (T/T c ) 2 , 
in which T c is to be understood as the mean field critical 
temperature of orbital ordering. We set T r — 150K in 
our calculation^. 

In the tetragonal phase, the orbital magnetic sus- 
ceptibility is found to be isotropic in the Fe-Fe plane 
and is almost temperature and doping independent for 
6.0 < n < 6.2(see Supplementary material C). This is 
reasonable since the orbital magnetic response is con- 
tributed by the whole band, rather than the electronic 
state near the Fermi level only. At 77 = 6.1, the bare or- 
bital magnetic susceptibility in the Fe-Fe plane is found 
to be about 7.3/ig/eV, which is enhanced to 10/Zg/eV 
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FIG. 2: Overlay of the band dispersion along the T — X and 
T — Y direction in the orthogonal phase. The orbital character 
is indicated by the color of the lines and the dispersion in 
the tetragonal phase is plotted in thin lines for reference. In 
the calculation we have set 772 = 30 meV. The dashed line 
indicates the Fermi level at n = 6.1. 



after RPA correction. This is already comparable to the 
observed total in-plane magnetic susceptibility at 200K 
in 122 systems, which is about 4.5 x 10 _4 erg/G 2 molAs 
(or 14/Xg/eV}^. As a comparison, the bare Pauli spin 
susceptibility is only about 2^|/eV. 

When a symmetry breaking perturbation of the 772- 
type is turned on, a two-fold modulation shows up in the 
in-plane orbital magnetic susceptibility. The angular de- 
pendence of the in-plane susceptibility at T = T c /5 is 
shown in Fig|3K. Here 4> denotes the angle between the 
X-axis and the direction in which the magnetic suscep- 
tibility is measured. The relative strength of the mod- 
ulation is about 2.6% before RPA correction and is en- 
hanced to 4.5% after RPA correction. The principle axes 
of the modulation are along the direction of the nearest 
Fe-Fe bond, which is just what we should expect from our 
model construction. The temperature dependence of the 
susceptibility in the principle axes are shown in FigJSjD. 
These predictions are in good agreement with the result 
of the recent torque magnetometry measurement^. Thus 
the magnetic anisotropy provide a realistic probe of the 
orbital ordering in the Iron-based superconductors. 

A robust prediction of our theory is the strong 
anisotropy between the in-plane and the out-of-plane or- 
bital magnetic susceptibility. At n = 6.1, the bare or- 
bital magnetic susceptibility in Z direction is found to 
be about 3.8/i|/eV, which is enhanced to 4.5/i|/el / af- 
ter RPA correction. This is only about the half of the 
value of the in-plane orbital magnetic susceptibility. We 
find the ratio between the in-plane and out-of-plane or- 
bital magnetic susceptibility is also almost temperature 
and doping independent for 6.0 < n < 6.2 and is always 
close to 2. According to experiments, both the in-plane 
and the out-of-plane magnetic susceptibility exhibit lin- 
ear temperature dependence with almost the same slope. 
However, the intercept of the in-plane magnetic suscepti- 
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bility is always much larger than that of the out-of-plane 
magnetic susceptibility 14-17 . This behavior can be eas- 
ily understood if we decompose the measured magnetic 
susceptibility into an isotropic component that is linearly 
temperature dependent and a temperature independent 
component that is anisotropic, or, 

X a (T) = xt+Xs(T). (6) 

It is then quite natural to associate the anisotropic com- 
ponent xt with the orbital magnetic response, which is 
essentially temperature independent. The isotropic com- 
ponent xs(T) should then be attributed to the spin mag- 
netic response, whose linear temperature dependence is 
still an unresolved issue in the fieldiS2r— 




4.(71) 




100 200 
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FIG. 3: (a)The in-plane modulation of the orbital magnetic 
susceptibility before and after RPA correction, (b) The tem- 
perature dependence of the RPA-corrected orbital magnetic 
susceptibility along the two principle axes of the orthogonal 
phase. 

In our calculation, we have used a five-band model de- 



rived from the band structure of the LaFeAsO system. 
However, the best known susceptibility data on single 
crystalline sample are all taken from the 122 system. It 
is thus better to perform the calculation with a material- 
spccihc band structure for the 122 systems. While this is 
an interesting possibility and should be pursued in the 
future, we note that the basic structure of the bands 
in both the 1111 and the 122 systems are quite simi- 
lar. Since the orbital magnetic response is contributed 
by the whole band rather than the electronic state near 
the Fermi level only, we expect the 1111 and 122 system 
to exhibit similar orbital magnetic response. Another 
way to improve our calculation is to use the matrix el- 
ement of L a calculated from first principle code, rather 
than approximating them with those in the basis spanned 
by the atomic orbital. However, since the form the ma- 
trix element is largely determined by symmetry, we do 
not expect such more advanced calculation to change the 
conclusion of this paper in a qualitative way. Indeed, we 
find that our results are not sensitive to the small varia- 
tion of the parameters 71,. ,5. 

In summary, we have shown that the orbital angular 
momentum of the conduction electrons in the Iron-based 
superconductors contributes significantly to the magnetic 
response of the system. In particular, the theory predicts 
that the orbital magnetic susceptibility accounts for more 
than 2/3 of the observed magnetic susceptibility at 200 
K in 122 systems. We show that the orbital magnetic 
response is sensitive to symmetry breaking in the orbital 
space, which makes it a useful probe of the orbital order- 
ing in these multi-orbital systems. A large and temper- 
ature independent anisotropy between the in-plane and 
the out-of plane susceptibility is predicted, which pro- 
vides a natural understanding on the behavior of the 
magnetic response of these systems. 
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and Tao Li is supported by NSFC Grant No. 10774187, 
No. 11034012 and National Basic Research Program of 
China No. 2010CB923004. We are grateful to K. Kuroki 
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I. SUPPLEMENTARY MATERIALS 

A. The form of the RPA kernel for orbital 
magnetic excitations 

The ten orbital magnetic excitation of the form O vv = 

<V, CT — c v t a Cv,o) are all time reversal odd and 
spin rotational invariant. In the absence of time rever- 
sal symmetry breaking they form a subspace within the 
space of all orbital excitations. It is thus sufficient to 
restrict our consideration in this subspace. 

The RPA correction to the orbital magnetic response 
is contributed by the inter-orbital Coulomb term, the 
Hund's rule coupling term and the pair hopping term. 
For example, the inter-orbital Coulomb term has the fol- 



lowing mean field decoupling(z/ > v), 

~ ^ ( C !,i/',<T C *,^< J ) C !,t>,o' C i,t / ',o' 

~l~ ^ ( c l,l/,o' c^ > ^y ^ o ')( c ll/',o■ c, ^ l/ . o ')• 

When expressed in terms of O v,v , we have 

(7 

Thus the RPA kernel is diagonal in the subspace of O u ' v . 
Following the same steps, it can be shown that the RPA 
correction contributed by the last two terms in Eq.Q 
cancels with each other. 



B. The form of the symmetry breaking 
perturbation in the orthogonal phase 

The form of the symmetry breaking perturbation in the 
orthogonal phase can be determined from the following 
group theoretical arguments. We first consider the form 
of the on-site symmetry breaking term. The point group 
around each Fe ion in the orthogonal phase is D 2 and has 
four one dimensional irreducible representations. Among 
the five MLWFs, |3Z 2 - R 2 ) and \XY) both belong to 
the identity representation, \X 2 — Y 2 ) belongs to the Bi 
representation, the linear combinations \XZ) + \YZ) and 
\XZ) — \YZ) belong to the B2 and B3 representation. 
Thus symmetry allowed on-site Fermion bilinear terms 
have the general form of 

+ ^5 ^(4,2,0- + C l3,<r)( C ^2, ( T + C 4j3 , CT ) 
i,(T 

+ PeJ2( c U,° ~ c l,3,<r) (ct.a.a ~ c 2 ,3, CT ) (7) 

In the tetragonal phase, the local symmetry around 
each Fc ion is promoted to D2d, which has four one di- 
mensional representations and a two dimensional repre- 
sentation. Among the five MLWFs, \3Z 2 — R 2 ) belongs to 
the identity representation, \XY) and \X 2 — Y 2 ) belong 
to the Bi and B2 representation, the linear combinations 
\XZ) + \YZ) and \XZ) - \YZ) form the two components 
of the two dimensional representation. For this reason, 
the bilinear form c^ ^c^i^, c\ A a c lA ,a, 4,5,a c h5,a, and 

cj 2 a Ci^,a + cj 3 a Ci^,a all belong to the identity repre- 
sentation of T>2d- When these symmetric perturbations 
are removed from Eq.©, we get the symmetric breaking 



7 



perturbation in the orthogonal phase, which now takes 
the form of 

AH = Al X] (4,2,^,3,0- + ct,3 )( T C i,2,<r) 
i,<j 

i,cr 

in which Ai = /3s — /?6, A2 = /S4. 

The above argument can be easily generalized to deter- 
mined the form the symmetry breaking perturbation on 
various bonds. In particular, we find there are in total 13 
independent symmetry breaking perturbations on near- 
est neighboring Fe-Fe bonds. The form of these terms 
are 

AH = AH S + AH P + AH d , 

in which 

AH S = ^^(ctj^Cf+^ + ctg^Cf+j^) 

2,<5,£T 
1,6,(7 

+ K 9 ^2(Vs4 ili<T Ci + s,2,o + P's4,l,a c i+S,3.<j) 
+ «10 ^2(Ps4,2,<r C i+S,5,<7 + Ps4^,a c t+S,5,a) 
+ Kll 22(P54,5,<T C i+6,2,<T + Ps4,5,a c *+S;3,a) 
+ K12 ^2 (l?s4,2,<T C i+S,4,'r - v's4,3,<j c i+&A,<?) 
+ «13 22(Ps4,4,<r c i+S,2,<T - v's4,i,<J C i+&,3,<7), 



If we restrict our consideration to the subspace 
spanned by the dxz and dyz orbital, then up to nearest 



and 



AH d = Kx ^ ds(4,2, a Ci+8,2,a + 4,3,cr C i+8,3,a) 
+ K 5 ^2 ^s4,l,<T C i+S,l,<r 



K6 ds4,4 . a C l+ SA, a 
K7^2d S 4,5, <T Ci+S,5,cr- 
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Here ps, p'5 are p-wave form factors, d^ is the d-wave 
form factor. The value of these form factors are illus- 
trated in FigH] 



FIG. 4: An illustration of the p-wave and d-wave form factor 
defined in the main text. 



neighboring hopping term, the only allowable symmetry 
breaking perturbation has the following form 

i,a 

+ 7/2 ^2 d s(4,2,a C i+S,2,° + 4,3,<T C i+W,<r) 
+ V3 ^2(4,2,* c i+S,3,a + 4,3,<7 C i+S,2,<r), 

in which 771 = Ai = /3 S - /3 6 , 772 = «i, r)3 = k 2 . 



C. The temperature and doping dependence of the 
anisotropy ratio 

Unlike the spin magnetic response, the orbital mag- 
netic response is contributed by both intra-band and 
inter-band process. As a result, the orbital magnetic re- 
sponse is much less sensitive to the variation of temper- 
ature and doping concentration of the system. In FigJSJ 
we present the temperature and doping dependence of 
the RPA-corrected in-plane orbital magnetic susceptibil- 
ity and the ratio between the in-plane and the out-of- 
plane orbital magnetic susceptibility. 

From the figure it is clear that both quantities have 
only small temperature and doping dependence. More 
specifically, the relative change of the in-plane orbital 
magnetic susceptibility for 6.0 < n < 6.2 is only about 5 
percent. The change in the anisotropy ratio is less than 
0.1. 



T(K) 



(b) 




T(K) 



FIG. 5: The temperature and doping dependence of the RPA- 
corrected in-plane orbital magnetic susceptibility(a) and the 
ratio between the in-plane and out-of-plane orbital magnetic 
susceptibility(b) for 6.0 < n < 6.2. 



